Generalized Universality of Higher Transverse Moments of Quark TMD Correlators 



M.G.A. Buffing/'El A. Mukherjee,2'3'[t] and P.J. Muldersi'l 



'Nikhef and Department of Physics and Astronomy, VU University Amsterdam, 
De Boelelaan 1081, NL-1081 HV Amsterdam, the Netherlands 
^Department of Physics, Indian Institute of Technology Bombay, Powai, Mumhai 400076, India 
'^Institut fiir Theoretische Physik, Umversitdt Tubingen, 72076 Tubingen, Germany 

(Dated: July 16, 2012) 



The color gauge-invariant transverse momentum dependent (TMD) quark correlators contain 
process dependent gauge links in the bilocal matrix elements. In this paper, we split these process 
dependent correlators into universal TMD correlators, which in turn can be parametrized in uni- 
versal TMD distribution functions. The process dependence is contained in gluonic pole factors, of 
which the value is determined by the gauge link. The operator structures of the universal TMD 
correlators are identified using transverse moments. In this paper, specific results for double trans- 
verse weighting of quark TMDs are given. In particular, we show that for a spin 1 /2 target one has 
three universal time-reversal even leading 'pretzelocity distributions', two of which involve double 
gluonic pole matrix elements and come with process dependent gluonic pole factors. We generalize 
the results for single and double weighting to TMD correlators of any specific rank, illustrating it 
for unpolarized, spin 1/2 and spin 1 targets. 

PACS numbers: 12.38.-t, 13.85.Ni, 13.85.Qk 



The transverse momentum dependent (TMD) correlators for quarks and gluons include not only the dependence on 
the longitudinal momentum fraction x but also the dependence on the transverse momentum of the partons (quarks 
or gluons). This enables one to incorporate spin- momentum correlations. The correlators in turn are parametrized in 
terms of parton distribution functions (PDF) and parton fragmentation functions (PFF) . The leading TMD distribu- 
tion and fragmentation functions in these correlators include besides the well-known spin-spin densities that survive in 
the coUinear case (where the transverse momentum is integrated) also spin-orbit densities. These provide for instance 
a natural interpretation for single spin asymmetries observed at high energies. The correlators are the nonperturbative 
objects that enter the description of high-energy scattering processes through a convolution with the perturbative 
hard scattering process. They constitute nonlocal matrix elements of the parton field operators. Collinear gluons 
exchanged between the soft and hard parts are resummed into the nonperturbative objects and show up as the Wilson 
lines or color gauge links that make the correlators gauge-invariant. For TMD correlators the nonlocality in the oper- 
ators is in the transverse direction as well as longitudinal (light-like) direction, and there is no unique way to connect 
the fields through the gauge link [H-g. The link depends on the process under consideration. In fact, in the case of 
single spin asymmetries it is the closing of the gauge link with the transverse gauge link at light-cone plus or minus 
infinity that plays a major role in distinguishing time-reversal even (T-even) and time-reversal odd (T-odd) TMD 
distribution and fragmentation functions . To study the p^-dependence it is convenient to look at the transverse 
moments, obtained by weighting the TMD functions with one or more powers of Pt 0, [!]• I* has been shown in 
Rcf. [2] that single weighted correlators relevant for different azimuthal asymmetries can be expressed in terms of two 
collinear correlators, the first one containing a T-even operator combination and the second one containing a T-odd 
combination. The latter involves a quark-quark- gluo n matrix element with vanishing gluon momentum and is known 
as the ETQS or gluonic pole matrix element 043- This matrix element appears in cross sections multiplied with 
a process dependent gluonic pole factor, which depends on the hard part of the process. Examples of such process 
dependent T-odd functions are the Sivers and Boer-Mulders function. The pr-weighted moment of the fragmentation 
correlators can also be divided into two parts similar to distributions, but here T-odd effects can come also from the 
fact that one has complex non-plane wave final states 3 • Tor fragmentation the gluonic pole matrix elements vanish 
and since for a given transverse moment there is only one specific operator combination, there is no process depen- 
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dence |15l - ll9l | . This is for example the situation for the Cohins fragmentation function. Generally, the dependence on 
the gauge link complicates the universality properties of TMDs as well as factorization issues. 

In terms of transverse moments, the study of TMD correlators becomes simpler. There remains a process depen- 
dence, but this is dealt with by process dependent gluonic pole factors, that depend on the hard part of the process. In 
other words, given a process one already knows which correlators are important. While single weighted moments are 
important for cos{ip) and sin((p) asymmetries, one needs higher p^-moments for cos{nip) or sin(n(p) asymmetries. For 
double weighted transverse moments, one looks at weighting with p"Pt- The pretzelocity TMD PDF hj^rp{x,p^) is an 
example for which double weighting is important. It has received some attention in the literature recently. It is a twist 
two chiral odd and T-even TMD distribution function. It contributes to the sin(3(/J~ ips) asymmetry in semi-inclusive 
deep inelastic scattering (SIDIS) f20','2l'l and to the cos{2(p+ipa — ipb) asymmetry in the Drell-Yan process involving two 
transversely polarized protons [22, 23]. In some models, such as the Bag model and spectator model, the pretzelocity 
distribution is shown to be related to the difference between the helicity distribution and the transversity distribution 
of the nucleon [2l|. This relation is not expected to hold in the presence of gluonic interactions. In this work, we 
analyze the double p^-moment of quark correlators taking into account the gauge link, and show that like the first 
moment these can also be separated into a T-even and a T-odd part. The T-even part contains three contributions, 
two of them coming from quark-quark-gluon-gluon matrix elements containing two zero momentum gluons, which are 
double gluonic pole matrix elements. The coefficients of these matrix elements depend on the gauge link U and are 
process dependent, showing that also the T-even pretzelocity PDF hf^\x,p'^) is nonuniversal. We will show that 
the pretzelocity function is a combination of universal functions, linked to the three possible T-even matrix elements. 
While these three functions themselves are by construction universal, it is a particular combination that appears in a 
given process with link dependence in the multiplicative coefficients. The appearance of three pretzelocity functions 
is a striking example of how the separation of the correlator into T-even and T-odd contributions is no longer enough 
to isolate the process dependent part of the correlator when higher transverse moments are involved. In addition 
to this, we will extend the transverse moment analysis to give definitions of universal p^-dependent functions of a 
definite rank. This will be done in general for targets with spin and illustrated for unpolarized, spin 1/2 and spin 1 
targets. 



II. FORMALISM 



A. Starting points 



The quark-quark TMD correlator is given by 



'^f^\x,pT;n) 



e'P-«(P|V',(0)%e]^.(O|P) 



(1) 



where we use the Sudakov decomposition p^ — xP^ + Pt + cn'^ for the momentum p^ of the produced quark. In this 
decomposition, is the momentum of the incoming hadron, which is in essence the leading light-like direction, while 
n is the conjugate light-like direction satisfying P • n = 0. The component a (x p ■ P along this direction is integrated 
over. The nonlocal matrix element ^[^'(SjPr; n) contains a process dependent gauge link C/[o.5]j connecting the two 
fields. The process dependence is in the path of the gauge link. For the two simplest possibilities, the and [—] 
gauge links, the gauge link runs from to ^ through plus or minus infinity along n, respectively. This is illustrated in 
Fig. [TJ More complicated gauge links can arise as well. We refer to Ref. [4] for a detailed description of these gauge 
links. After integration over transverse momenta, one has the quark-quark collinear correlator. 
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FIG. 1: The gauge links [+] (a) and [— ] (b) running from to ^ with £ ■ 
transverse separation are nonzero. 
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where the gauge link is reduced to a straight-hne gauge link or Wilson line, which runs from to ^ along n. Since the 
quark-quark correlators cannot be calculated directly, it is common to make a parametrization that contains TMD 
or coUinear PDFs, respectively. In the TMD case, there are for a spin 1/2 nucleon eight leading contributing terms 
in the parametrization of the TMD correlator , 

+ hfj{x,pl) 75 $r + /.^-fl ^ + ^hi^''\x,pl) |}|, (3) 

with the spin vector parametrized as S"^ = S^P^ + S!^ + Sj^nf^ and shorthand notations for g^^' and hf^'^\ 

9u\x^Pt) = Sr.gfi]{x,pl) - ^'^ gfrix^pl). (4) 

The TMD distribution functions in this parametrization depend on x and p^ — —p^ — — |prp- The leading contri- 
butions in the correlator all have a f factor and are distinguished by different azimuthal behavior for the transverse 
vectors such as pr and St- The correlators and the TMD distribution functions in the parametrization also depend 
on the gauge link. Time-reversal relates the functions in ^^^^ to those in $1^ \ where C/* is the time- reversed gauge 
link, which means interchanging the running via light-cone plus or minus infinity. For the functions and hj^ one 

has f^rjy^ — —fiT^ \ a property that is referred to as naive T-odd. Each TMD has either zero, one or two factors 
of Pt as a pre-factor. This will play a role when integrations over transverse momenta are considered. It is actually 
useful to use in the parametrization irreducible (symmetric and traceless) tensors in the transverse space, 

Pt, Pt — Pt Pt I^Pt yT ■ [O ) 

Just integrating (without weights) Eq. [3] over transverse momenta, only the contributions without pre-factor of Pt or 
traces survive, yielding 



Hx) = |/i(x) + S,g,ix)-f5 + h,{x) 75 



(6) 



at the leading twist two level. Here gi{x) is the integrated version of gf^{x,PT) and hi{x) is the p^^-integrated version 
of H^^^XjPt) = h^^{x,PT) + hf^^^-^\x,PT) including a trace term, which involves functions weighted with powers of 



pI/2M^ = pII2M^, in general 



/ 2 \ " 

/.^."^(^,P^)= (^) f...{x,pl). (7) 



The integrated functions /'"-'(a;) are usually referred to as transverse moments, but we will extend this name to 
azimuthally averaged functions that still depend on p^. The collinear PDFs in Eq. [6] are independent of the gauge 
link U . In other words, all operator definitions of these collinear PDFs have a unique straight-line gauge link. 

The behavior of (TMD) PDFs under time-reversal can be studied. The functions and are time-reversal odd 
(T-odd), while the remaining six functions are time- reversal even (T-even). Similarly, one can look at the behavior 
of the matrix element(s) under time-reversal. Using the fact that the simplest gauge links for quark correlators, the 
and [— ] gauge links, are a time- reversal couple, one can construct T-even and T-odd TMD correlators 

<i>(T-<=n)(a;,p^) = i($[+l(x,p,)-H$H(:c,p^)) , (8a) 
^'^'^-°'"'\x,Pt)^\{^^+Hx,Pt)-^^-\x,Pt)) . (8b) 

For the unweighted integrated case the separation between T-even and T-odd objects would be trivial, since the 
and [— ] gauge links are identical after integration over transverse momentum. As a result, $(a;) = (x) and 

^(T-odd)j'2.^ _ Pqj. ^jj^g transverse momentum weighted case both functions are important. One thus is tempted to 
identify the TMD functions and to the T-odd correlator and the other TMD functions to the T-even correlator, 
in which the T-odd ones acquire process dependence. The situation will turn out to be more complex, which is most 
easily demonstrated by looking at transverse momentum weighting. 



4 




X-Xi 



X— Xl 



(a) 



(b) 



(c) 



FIG. 2: The correlators (a) $^1(2; ~ xi,xi\x), (b) ^a{x\xi, x — xi) and (c) "I>A(a;, —xi \x ~ xi). Note that the diagrams in (b) 
and (c) are equal to each other. 



B. Single transverse weighting 



In the py-weighted case multiple matrix elements appear, since the transverse weighting gives rise to a derivative 
that not only acts on the fields, but on the gauge links as well. Weighting with can be rewritten in terms of 
two contributions, which upon pr-integration only depend on x and depend on the link just through a gluonic pole 
factor 0, 



nix)+7:d^^mx 



(9) 



The matrix element <&Q(a;) is referred to as the gluonic pole matrix element. All matrix elements are built from multi- 
parton twist three operator combinations, illustrated in Fig. [2l The relevant ones in the final result of a calculation 
involve gauge- invariant operators iD" and F"" rather than A", 



'^Dtji^ - xi,xi\x) 
^Fiji^ ~ xi,xi\x) 



d£,-Pd'n-P 

(2^)2 
d^-Pdrj-P 
(27r)2 



,ipi-V+iip-Pl)< 



{P\ij,{0)U[o,r,]^D'^{v)Ui,^^]MO\P) 
{P\^^m[o,n]Friri)U[n,iMO\P) 



LC 



LC 



(10) 

(11) 



The matrix elements showing up in Eq. ^ are related to these multi-parton correlators, to be precise we need the 
matrix elements 



^d{x) = J dxi ^"^yix - xi,xi\x), 

= I dxiPV—^T'(x~xi,xi\x), 

J Xl 



(12) 
(13) 
(14) 



By using the Eqs [8al and I8b| one finds that <&g(a;) is a T-even matrix element, involving the and the latter 
being the principal value integration over a correlator involving the gluon fields defined in a T-invariant way. The 
gluonic pole matrix element $q (x) , in order to distinct it from $ p indicated with an index G, is T-odd. 

The choice of notation for the arguments of multiparton correlators, where the produced quark and gluon on the 
left side of the cut have momentum fractions x — Xi and Xi and the incoming quark on the right side of the cut has 
a momentum fraction x, is illustrated in Fig. [2^. Despite the fact that the momentum assignment in the correlators 
in Eqs llOlfni are overdetermined, it has the advantage that it is more transparent in our forthcoming generalization 
to higher weightings. Furthermore, because of the absence of T-ordcring, one can move a gluon through the cut by 
changing the sign of the momentum. Under hermiticity one finds that the correlators in Eq. [14] have the behavior 
7o$^(x - a;i,xi|x)7o = ^a{x\xi,x - xi). 

The weighting with transverse momenta can also be analyzed by studying the parametrization in PDFs. For single 
Pr-weighting, only PDFs with one pre-factor of pr in the parametrization in Eq. |3] survive. T-even PDFs contribute 
to the ^g{x) matrix element, while T-odd PDFs contribute to the ^q{x) matrix element, see Ref. for a detailed 
study of this. Since the T-odd matrix element comes with a process dependent pre-factor, it can be seen that for single 
Pr-weighting, the behavior under time-reversal can be used to identify the process dependent parts. In literature, it 
has become common to use notations like or f^j^^^^ for these functions. This suggests that there are many 

different versions of specific functions, which are obviously not universal, one for each process with a different gauge 
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TABLE I: The values of the gluonic pole pre-factors for some gauge links needed in the pr-weighted cases. 
Note that the value of Cq^' is the same for single and double transverse weighting. 



link. In fact there is a universal transverse moment relating all link dependent ones 

ft^'^^^\.)=cffti'\.). (15) 

Although the only difference for the single weighted case is just the numerical pre-factor that for simple processes 
is just +1 or —1, we will show in the next section that for the double weighted case the situation becomes more 
complicated and one actually gains a lot by this different notation. But even for single weighting there is a clear 
advantage using Eq. [TSl because it states that there is a universal function with calculable process (link) dependent 
numbers rather than an infinite number of somehow related functions. For some gauge links, these numbers are shown 
in Table n Here C/Pl is the Wilson loop J/I-l ' J/W. 



C. Double transverse weighting 



In order to evaluate the double transverse weighting we need to consider matrix elements like 

J ZTT ZTT ZTT 

X (p, s\m ul;l]Fr{v')ul:;l,^F-^{v)ul;]^^ m\p, s) 



(16) 



LC 



among others, where LC indicates that all transverse components and rt-components of the coordinates are zero. 
Besides this matrix element one needs $_d_f, ^fd and ^dd as well as bilocal matrix elements, obtained by direct 
or principal value integrations over these matrix elements (as in the case of single transverse momentum weighting) 
or gluonic pole matrix elements, where xi or X2 or both are zero. Explicitly, the matrix elements are discussed in 
Appendix [X] 

The actual weighting of the gauge link dependent TMD correlator ^^^^{x,Pt) gives 



lu] 



(x) ^ J dVpi%.'^$[^l(^,P^) 



^ifHx) + .4^1 ($r/>(x) + %-^Hx)) + .^Cgl, ^l?^lix) + .^Cl^, (,). (17) 



For the correlators containing two (or more) gluon fields like the one in Eq. I16[ one must distinguish the different 
color structures for the correlator, hence a summation over the color structures c. For double weighting, there are in 
the double gluonic pole part two possible color structures related to the appearance of the color traced Wilson loop 
Ttc{U^^^). The differences between the two different correlators ^c'Q\ix) are made explicit in Appendix El Just 
as for the single weighted case in Eq. [HI the structures with one or more partial derivatives denote differences 
between correlators with a covariant derivative minus a correlator with a principal value integration, e.g. $^^'^^(2;) — 

^o'c^i^) ~^M3^ (■'')■ ^'-'^ completeness, they are given in Appendix lAl Since the weighting is done with the symmetric 
combination, we have symmetrized in the indices, which should not infiuence the result. We also omitted the Dirac 
indices on the fields. The precise form of all correlators in terms of matrix elements can be found in Appendix [XI 
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The only leading twist TMD PDF that contributes is the Pretzelocity hf^^XjP^). Since this function is T-even, 

its (double) transverse moment could be associated with both the matrix elements ^Qif\x) and 7r^$g^|^(a;). The 
giuonic pole matrix elements come with gauge link dependent pre-factors, so the Pretzelocity function as it has been 
defined in literature up to now is not universal. The double giuonic pole factor gives the gauge link dependence and 
one must identify the gauge link dependent function as the sum of three functions, 

htr^Hx) = ht^'^^-\x) + + C!^k.htr'''\-)^ (18) 

where the functions hf^'^''^'^\x), hf^'^^^'^^\x) and hf^^''''^^\x) are universal. For simple processes like SIDIS and 
Drell-Yan with Cgg,i = 1 and Cgg,2 = one finds just the sum of two Pretzelocity functions. For processes with 
a more complicated gauge link structure, other combinations involving three functions will appear, as can be seen 
m Table H The double -weighted results also show that for higher transverse moments, and hence also for the full 
PT-dependent treatment, a separation in T-odd and T-even functions is no longer sufficient to isolate the process 
dependent parts. 



III. DEFINING TMDS 



For the definition of a TMD correlator parametrized in terms of PDFs (or PFFs) depending on the collinear 
fraction x and transverse momentum it is important to keep in mind the role of x and Pt- These are identified 
with kinematic variables in a high-energy scattering process. This is most well-known for x, which in a SIDIS process 
is identified with the Bjorken scaling variable. In the same way the transverse momentum can be identified, e.g. from 
the noncollinearity of produced hadrons or from jet-jet asymmetries, even if for transverse momenta the identification 
is usually contained in a folding of transverse momenta of several hadron correlators. So for purposes of further 
analyzing we assume that we know that a hadron correlator depends on 'measurable' x and pr- After integrating out 
the 'unmeasurable' transverse momenta of the soft gluons contributing to the Wilson lines, this is at most a single 
transverse momentum. To incorporate this we have made sure that in the definition of multi-parton correlators these 
components correspond to the total momentum p flowing between a (soft) hadron correlator and the hard process. 

The leading relevant TMD operator structure for our considerations thus will be of the generic bilocal form 



(27r)3 



(19) 



C-n=0 



where U is one of the possible gauge links for TMD correlators. We can define in this way ^g{x,pT) = ^%{x,Pt) 
^'X{x,Pt), as well as correlators $?fn, $ 



Q/3 



where one in particular for ^gg must take care 
of the color structure. 

Our identification of operator structures and TMD functions in the parametrization of correlators depends on the 
comparison of moments in x and pr, even if such moments in real life are limited by kinematics of the process. This is 
well-known, but nontrivial, for the collinear dependence, where the moments can be related to local matrix elements of 
quark and gluon fields. All these operators have the same twist (canonical dimension minus rank of Lorentz indices), 
which means they contribute at the same order of the hard scale. The x'^~^ Mellin moments correspond to expectation 
values of leading twist operators of rank TV. These local matrix elements have a calculable scale dependence governed 
by the anomalous dimension of the local operator. The scale is usually identified with the kinematic limit such as 
the exchanged momentum in SIDIS. The a;-dependent functions can be reconstructed from the Mellin moments. 
Their scale dependence then is obtained by folding them with splitting functions, of which the Mellin moments are 
precisely the anomalous dimensions. 



A. TMDs of definite ranlc 



For the p^-dependent functions we follow a similar procedure. An expansion of TMDs involves the symmetric 
traceless tensors p"^ '"'" of rank m. These traceless tensors satisfy 

dVK"-"'"PTn...^™/...(^,p') « /.^."'W- (20) 
Actually, it is sufficient and for our purposes desirable to integrate only the azimuthal part, 

^P^'■■■"-PT^,...^^f...{^,pl) OC I^:"'\x,pI). (21) 
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The r.h.s. of these equations contain the transverse moments defined in Eq.[7|as well as constant tensors without 
azimuthal dependence. Writing the following parametrization, 

V [x,Pt) — ^[x,p^)+ttUq 'i'd^^PT) + ^ ^gg,c jyp ^ggA^'Pt) + ^ '-^ggg,c j^p ^ggg,c\^^Pt) + ■ ■ ■ 

^ G Jp ^{dGj^^^PT) + ^ <-'GG,c ^{OGGjA^'PT) + • ■ • 

+ '-^^^aai-yr) + .-., (22) 

we reproduce the moments. Note that in Eg. 1221 depending on the gauge link there are multiple contributing color 
structures for terms with two or more gluonic pole terms, hence the inclusion of the summation over these color struc- 
tures c. For the collinear correlators ^^qq ^{x) this is discussed in Appendix [XI for the TMD correlator ^{x,p'^) it 
is discussed in Appendix[Bj The operator structures on the r.h.s. thus are the ones appearing in an angular expansion, 
in which the azimuthal dependence is made explicit. The combinations like ^{og} indicate symmetrized combinations 

^{dG} — ^dG + ^Gd- Upon (p-integration only the structure ^(x,p^) survives. Hence, we identify this as the rank 
zero TMD correlator. 



<^ix,p^) = |/i(x,p^) + Sr^gi{x,p^)-/5 + hi{x,p^)-f5 PtJ^, (23) 

where the TMD correlator hi{x,p^) rather than hiT{x,p'^) appears (see also the remark following Eq. [6]). Next, we 
look at the weighted expressions before pr-integration in order to explicitly identify further TMD functions, 

g$[^l(.,p.) = ^^x,pl)~^'^,^'\x,pl)-nd^^<^%^'\x,pl) 

+ E^^l{:,,pl) + nCl^^^<i>hix,pl) + ... . (24) 

The first term is obviously a term that needs to be there because we already identified a nonzero rank TMD 
correlator. The next two terms involve the T-even and T-odd transverse moments of <&g(a;,p^) and $g(x,p^), 
respectively, which survive (/?- integration. The other terms in Eq. [23] contain higher rank tensors in pj,. Comparing 
the unintegrated expression with the parametrization for a spin 1/2 target one thus immediately identifies in addition 
to the rank correlator in Eq . [23] two rank one TMD correlators, 

— $a(a;,p^) ^ <hj^Lix,p^)SL^^ ~ 9it{x,p^) ^^-f5>^, (25) 



P 



Tl 



^hi^,pl) = { - fM-,pl) '-^P^ + ^hiix,pl) 



the first one being T-even, the second one T-odd. 

Before going to double weighting, it is useful to realize that in Eq. [22] one does not need to subtract trace terms 
from the second rank correlators. This is automatic because of the use of the tensor p"^ . We can write 

Jp'^ldG} ~ [-Jp- - ^9T^J j ${aG} -^^'^{dG} + '^{0-G} " [^{dG}-2 9T *{9G} J ■ (27) 

In Eq.[27]we introduced the notation d-G in the subscript of one of the correlators to indicate that these two operators 
in this correlator have been contracted. Also for double p. -weighting we write down (selected contributions in) the 
unintegrated result starting with Eq. [51] We find for the double weighted result, 

^ (^tPt) — j^p 'i'(x,PT) ]\^ d \^->Pt) 2^ G 9t ^ G \^^Pt> 

- ^p'^ ~^T\-^pI) - ^T\-^pI) + • • ■ 

- I af {^^9%ix,pl) + Trdc^ ^fLyi^.pl) + vr^^G.c *G.G,c(^'P' )) + ■ • ■ ' (28) 
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an equation that after symmetrizing in a and /3 can be used to identify the remaining TMDs for a spin 1/2 target. We 



have omitted terms with rank 1 tensors in multiplying ^qqq, etc., as well as terms with rank 3 or rank 4 tensors 



like and p^^^j ■ The last line in Eq. [28l containing terms that arise by contraction of the indices i and j in the 

explicit rewriting of the product p'^PrPrij through the contribution gr^gnj, is included for completeness. Note that 
taking the trace —hgrap gives 



p 



2 

1 



_ „2 N _ l^clU] $(2) „2^_ 1 2^[C/] ^(2) / 2 N , 



where the terms between brackets come from the correlators with contracted operators in Eq. [28l This shows that 



- ^ {pi" ^T\^y.) - trace) - 2^ TTClfl [pi- <^§^'\x,pI) - trace) + . . . 
+ {^fe^'\x,pl) - traces) + ^cg'l - traces) 

+ T^'^G.c ('J'gg.c ) - traces,) + . . . , (30) 

illustrating how the projection with the properly symmetrized traceless second rank tensor p%^ gives the properly 
symmetrized traceless TMD structures. The terms without azimuthal dependence are identified with the rank two 
TMD correlators, which for a spin 1/2 target are parametrized as 

ddy^^Pr) — ''-IT [^iPtJ ^2 2' 

]y,j2^GG,l\^^PT) — [X,Pt) 2' ^ ' 

^<^' fx «2^ :Lh^^B2), 2.PT^2S^J5l^ f_ , . 

j^2^GG,2\^^Pt) — ^2"'IT [X,Pt) 2' ^ ' 

J^noG}i-^Pl) = ^- (34) 

The last TMD correlator in these equations is a T-odd rank 2 TMD correlator, which is not present for a spin 1/2 
target. The result can also be summarized as the existence of three universal Pretzelocity functions hf^'^\ hf^^^^ 
and /i^^^^' and a gauge link dependence given by 

f^lT H^iPt) — ^IT \^iP't) ^GG.l^lT \^tPt) ^GG,2^1T \^iPt)- (35) 

This shows e.g. that hf^^^^^\x,pl,) — hf^"^^ (x,p^), but that for other processes (with more complicated gauge links) 
other combinations of the three possible Pretzelocity functions occur. In asymmetries involving p"'^-moments of the 
quark TMD correlator contributions from all four correlators can appear. In particular we find for a transversely 
polarized spin 1/2 target three Pretzelocity functions, as was already established in Eq. [TH] For a spin 1/2 target our 
treatment is complete, since there are no higher rank TMD correlators such as Pnjk^ ^gg{x , p^) ■ In the case of a spin 

1/2 target the Pretzelocity TMD functions /i^^^^^ actually was referred to as junk TMD in Ref. Q. 

We want to summarize our results in this section in tabular form. We first represent the contributions in Eq. [52] in 
Table nil The assignment of the TMD PDFs for an unpolarized and polarized spin 1/2 target has been discussed in 
this section and is summarized in Tables IIIII and IIVI For the corresponding fragmentation functions the assignments 
are different, since gluonic pole matrix elements vanish in that case and all functions are assigned to ^d...d 

operator structures. The assignments thus are as in Tables |V] and IVII The T-odd TMD PFFs (such as the Collins 
function H^) are due to the fact that the definitions of fragmentation functions involve non- plane wave states or 
equivalently a hadronic number operator, which are not invariant under time-reversal. There, thus, is only a single 
(T-even) function H^{z, fcj) appearing in the parametrization of the correlator A'^g{x,p'^). 
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GLUONIC POLE RANK 





1 


2 


3 








'^^C'ggg.c '^GGG,c 






^'C^GG,c3>{aGG},c 












'I'aee 









TABLE IL The contributions in the TMD correlator for correlators ordered in columns according to the 
number of gluonic poles (G) and ordered in rows according to the number of contributing partial derivatives 
{d = D — A). The rank of these operators is equal to the sum of these numbers. Their twist is equal to the 
rank + 2. 



PDFs FOR SPIN HADRONS 


/i 


hi 

















TABLE IIL The assignment of TMD PDFs 
for a spin or unpolarized target to the quark 
correlators as given in Table |II] involve at most 
rank 1 TMD correlators. _ There is no T-even 
function corresponding to $g. 



PDFs FOR SPIN 1/2 HADRONS 


ffi, hi 


flT 










"it 







TABLE IV: The assignment of TMD PDFs for 
a polarized spin 1/2 target to the quark corre- 
lators as given in Table |TI] involve at most rank 
1 TMD correlators for longitudinal polarization, 
while they involve also rank 2 TMD correlators 
for a transversely polarized spin 1/2 target. 



PFFs FOR SPIN HADRONS 


Di 



















TABLE V: The operator structure of quark 
TMD PFFs for spin or unpolarized hadrons. 
All gluonic pole matrix elements vanish. 



PFFs FOR SPIN 1/2 HADRONS 


Gi, Hi 






GiT, -ffiL, -Drr 






Hit 







TABLE VI: The operator structure of quark 
TMD PFFs for polarized spin 1/2 hadrons. Glu- 
onic pole matrix elements vanish. 



B. Results for spin 1 hadrons 

Extension to higher spin targets is straightforward. We illustrate this by giving in Table IVIII the assignments for 
spin 1 tensor polarized TMD functions. These were first given in Ref. [25j . The (slightly updated) parametrization 
of the TMD correlator for the TMD PDFs for a tensor polarized target are given in Appendix [C] as well as the 
parametrization of the TMD PFFs in Appendix ID] From these tensor polarized spin 1 contributions, the /^yj-' (-^'Pr) 
and hf^{x,p'^) can be written as a combination of multiple universal PDFs, multiplied with process dependent 
gluonic pole factors, 

/iTT'('^'Pr) — flTT^ (■'^tPt) ^ ^GG,c flTT \^tPt)j (36) 

Note that the (xjp^) is a rank 3 object, for which all contributing universal functions are multiplied with a 
process dependent pre-factor. A special case is the T-odd TMD PDF Hilt, which is forbidden because of time- reversal 
invariance. Following Ref. [25||, this rank TMD PDF is defined as the combination hiLT{x,p^) = h'-^i^rp{x,p^) + 

^iLT (^tPt) is shown as the wiped-out function in Table IVlIl It shows a nice feature of our TMD functions of 
definite rank. In the first column only T-even TMD PDFs are allowed, in the second column only T-odd ones, etc. The 
first victim of the application of time reversal invariance for leading quark TMDs, thus, is hiLTix,p'^), a (T-forbidden) 
transversely polarized quark distribution function in a tensor polarized hadron. Note that the rank 2, T-odd function 
^ilt(^'Pt) allowed. The only rank function for a tensor polarized spin 1 target thus is fiLL{x,p'^), introduced 
as the distribution bi in Ref. [26]. 

For fragmentation functions, gluonic pole contributions all vanish and only the first column survives. The 
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PDFs FOR TENSOR POLARIZED SPIN 1 HADRONS 


flLL, ^bfffpi 


hlLL, giLT, hiTT 


ABc) 
J ITT 




flLT 


hlLT^ 9lTT 






AA) 
J ITT 


1 TT 















TABLE VII: The operator assignments of TMD PDFs for a tensor polarized spin 1 target require operator 
structures up to rank 3. There are several different functions /itt(2;,Pt) and h'^j,j,{x,p%). 



PFFs FOR TENSOR POLARIZED SPIN 1 HADRONS 


Dill, Hilt 








DiLT, Hill, Gilt, H\tt 








DiTT, H^lt, Gitt 






Hvtt 





TABLE VIII: The operator structure of TMD PFFs for a tensor polarized spin 1 target requires operator 
structures up to rank 3. 



parametrization of the higher rank correlators contain the T-even and T-odd TMD fragmentation functions. The 
fragmentation functions describing fragmentation into a tensor polarized hadron are given in Table IVIIII 



C. Bessel weights 

We note that the TMDs f^"'\x,pl) of a given rank do not contain operators of definite twist. This is only 
true for transverse moments /'-'"■'(a;) after pr-integration. The TMD correlators of definite rank appearing in the 
parametrization in Eq. [22] only arc integrated over azimuthal directions. The rank just refers to the azimuthal 
dependence of the correlators in the full correlator <^^^\x,Pt)- 

Using that for a given rank m, there are two independent combinations p^i ' *™ oc |pr|™ exp(±im(pp), it is equivalent 
to consider 



VtII . . .i-n 



<i>'\-'"^ix,pl) or $(™/2)(2.^p2)g»mVp^ (38) 



where — (— p^/2A'f^)"'/^ $...(a;,p^) assures the appropriate small pr-behavior. A suitable normalization 

of the correlator has to assure that $'™-'(a:,p^) reproduces the collinear transverse moments upon integration, 

/>oo 

2TT\p^\d\pr\ ^^"'\x,pI). (39) 



Knowing the correlators in Eq. [35] to be Fourier transforms of nonlocal matrix elements in transverse space, it is 
natural to write the appropriately weighted TMD PDF in their parametrization as a Bessel transform, 

/>oo 

/ db./\^bJ^{\pr\b)f^:"^/'Hx,b), (40) 



such that Z^™/^^ (a;, 6) exp(im(pf,) is the (two-dimensional) Fourier transform of /'"'^^•'(a;, |pr|) exp{irrnpp). Bessel 
weightings are extensively studied in Ref. ^7| . 



IV. CONCLUSIONS 



In Eq. [12] we have presented a parametrization for TMD quark correlators that distinguishes different azimuthal 
dependences by writing an expansion in terms of irreducible tensors in the transverse momentum multiplied with 
correlators depending on x and p'^. These correlators contain tensors describing the polarization of the target and 
TMD functions depending on x and p^. The rank of the irreducible tensors in transverse momentum space also defines 
the rank of the correlators and TMD functions multiplying this tensor. The field theoretical expression for the quark 
correlator basically has two quark fields connected by a gauge link. The operator structure of the TMD correlators of a 
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definite rank contain Dirac gamma matrices, derivatives or gauge fields with transverse indices in color gauge-invariant 
combinations and a definite rank. They are structured in a similar way as higher twist operators in the collinear case. 
Each independent operator combination defines a particular TMD function. For leading twist operators the relevant 
'transverse' operators are either gluonic pole (G) or partial derivative {d = D — A) or combinations thereof. A special 
feature of these operator combinations of definite rank is that they are either T-even or T-odd after extraction of a 
gluonic pole factor without having to perform weighting and integration over transverse momentum. 

Using the parametrization of TMD correlators in Eq. 1221 one finds that rank and rank 1 contributions are similar 
to previously used definitions for T-even and T-odd contributions, such as e.g. obtained by combining 'opposite' gauge 
links in Eq. [H Just as the collinear transverse moments, these TMD functions are universal functions, multiplied 
with a process dependent pre-factor, rather than nonuniversal gauge link dependent functions. By explicitly and 
systematically looking at all contributing functions, one finds for an unpolarized target two TMD quark correlators, 
the first being of rank 0, containing the TMD PDF /i(x,p^). Looking at the operator structure, it is interpreted 
as the momentum distribution of quarks. The second unpolarized TMD correlator is a T-odd gluonic pole matrix 
element of rank 1. It contains the Boer-Mulders TMD PDF hi{x,p'^). 

For a polarized spin 1/2 target there are an additional eight TMD correlators, containing rank 0, rank 1 and rank 2 
contributions. The two rank correlators contain the TMD PDFs that are interpreted as the well-known polarized spin 
distribution functions in longitudinally or transversely polarized targets. For a longitudinally polarized target, there is 
a T-even rank 1 TMD correlator <i>g(a;,p^) containing the worm gear function h-^^{x,p'^). For a transversely polarized 
target, there exist one T-even TMD correlator ^q{x,p^) and a T-odd rank 1 TMD correlator ^c{x,p^) containing 
the worm gear function giT{x,p'^) and the Sivers function f^{x,p'^), respectively. The three TMD correlators of 
rank 2 appear in the T-even correlator <I>gg(x,pJ) and the two T-even double gluonic pole correlators ^Gc dx^p^), 

giving rise to the Pretzelocity functions hf'^'^\x,p'^), hf^^'^\x,p'^) and hf^^'^\x,p'^). These functions in general 
both show up in particular azimuthal asymmetries but with gauge link dependent pre-factors, where the gauge link 
in turn depends on the process. As for the functions themselves and in particular their interpretation, the function 
/i^^"*'* is related to the quark structure of a nucleon, while the functions hf^^'^^ are the ones involving quark-gluon 
correlations. 

For a spin 1 target one finds apart from the above mentioned TMD correlators additional correlators because one 
also has the possibility of tensor polarization. The full list of TMDs has been given in Table I VIII including rank 3 
contributions coming with process dependent pre-factors. Rank 3 contributions like this are specific for targets with 
spin 1 or higher. 

The procedure for defining universal TMD correlators of definite rank can be extended to gluon TMDs and to 
higher twist situations. The extension to gluon TMDs will be presented in a forthcoming publication. The situation 
for higher twist TMDs is complicated by the fact that the lowest twist operators that contribute to the TMDs not 
only contain two quark fields or two gluon fields, but they also contain additional (gluon) operators with transverse 
directions, D" and F"", no longer in the combination 9". 

As a final advantage of the universal TMD correlators we mention that, although their nonlocal operator structure 
is of the form in Eq. [19] with a particular gauge link U, the TMD correlators of definite rank have definite T-behavior 
(even or odd) and are independent of the gauge link. The ?7-dependence is in the gluonic pole factors and the color 
structure of the operator combination. Thus one can study the universal TMD correlators, for instance in lattice 
calculations, by using just the sum and difference of the simplest J/I+l and U^^^ staple links. 
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Appendix A: Double weighted collinear matrix elements 

Double pr-weighting is worked out in the same way as the single weighting. In evaluating the weighting one can 
actually choose the derivatives to work on ^ or or mix these, but that does not matter for the final answer. Since 
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we just are interested in the general structure, we work with the multi-parton element with all gluons on the left side 
of the cut. For the gauge link U^q^^^ we note that (writing only the relevant pieces of the matrix elements) [sO] 

^"t^[0,(_] - "^[0,00] ^"T^[OT,iT] ^[oa,(.] - ^[0,oo] ^ [Ot '^[oo,^]' \^^) 

which is then further evaluated using 

iD'i Ul^^^^ . . . m = f/g^g [iD-iO ~ A%{i) + ^G""(0) . .-m- (A2) 
The A%{i) and 7rG""(^) are defined as 



7rG""(0 = - / di^-P C/f'^";]G""(ry)C/[';;^,, (A4) 

^ J — oo 

with e(C) being the sign function taking the values 1, —1 and 0. Note that G""(^) = G""(^-P, ^r) does not depend on 
^•n. The above described method for calculating the structure of the matrix elements for a single transverse weighting 
can be extended to higher transverse weightings by repeated application of Eg. IA2I For example, for the C/^"*"! gauge 
link this implies 

= ^[iei - ^?(o) - • ■ • m + t/f::!,c] (^G""(c) ^G"^(c)) . ..m 

(A5) 

After integration over the resulting correlators can be rewritten in terms of color gauge-invariant multi-parton 
correlators as was done for the single weighting. In this case one needs the correlators of the form 

x(F,5|^(0)C/f„"l,,,O^(770<',]OL(^)<'a'^(OI^,^)jLC, (A6) 

with Ofy and Ojr hermitian operators like il?" and/or F"". In general more than one color structure is possible. 
As an example, for the operator combination tpFFip, one can have the two distinct color configurations 

c = 1 : Trc [FF^Pi^] = ^FF^P = V'''^''"'^"'"'/'", (A7) 
c = 2 : Trc [FF] Tr^ [tip] = V'V' Trc [FF] = i^"" i;'' F"''' F"' " . (A8) 
An example of a correlator one needs is 

<d(^) = |^e-(«-^\P,5|V^(0)C/fo"l,zZ??(0*i5^(0^(OI^,^)jLC 

= J dxidx2 ^%^^{x - xi - X2,xi,X2\x). (A9) 

Others correlators involving A" (^) or vrG"" (^) , given in Eqs IA3I and IA4[ are given by 

<t>'f^{x)^ I dxiPY— I dx2PY—^pp{x - XI - X2,xi,X2\x), (AlO) 

= / dxidx2Py—^p%{x-xi—X2,xi,X2\x), (All) 

^'^^a{x)= [ dxidX2PY—^lU'^-Xi-X2,Xi,X2\x), (Al2) 
J X2 

I dx2<i>fD{x-X2,0,X2\x), (A13) 

*Sg(^)= f dxi<i>fp{x-xuxi,0\x), (A14) 
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<^GAi^)= [ dX2PV—<^pUx~X2,0,X2\x), 

*AG(a^)= ! dxiVV—<^fp{x-xi,xuQ\x), 



(A15) 

(A16) 
(A17) 



Using for the correlators, just as for the single weighted case, the notation for the correlators with covariant 
derivative minus a correlator with a principal value integration {iDt — At), implies 



















-^"g'a 



The second transverse moment in terms of the coUinear functions then is (symmetrizing in a and /3) 



(A18) 
(A19) 
(A20) 

(A21) 



which is the result given in Eq. [17] with C^^' — 1, c\^q i — ^ ^-nd 2 = 0. For other gauge link structures, similar 
calculations can be performed. 



Appendix B: Double weighted TMD matrix elements 

The leading relevant TMD operator structure for our considerations referred to in Eq. [19] is bilocal. 



(27r)3 



(P|^(0)C/[o,e]0(e)|P) 



(Bl) 



The nonlocality, however, involves a transverse separation, hence the gauge link t/[o,^] in general can be complicated. 
For the two cases c = 1, 2 one now finds for the gauge link C/j^'^j'''' the nonlocal structures 



c= 1 
c = 2 



Tr 
Tr 



[/f°l]G(OG(e)] Tr, [C/W, mm] - ^(0) t/fotJ] ^(^Tr, [u]^\f{i)G{i)U\^l^ 
Appendix C: Parametrization of the spin 1 distribution correlator 



(B2) 
(B3) 



The parametrization of a distribution correlator for a spin 1 hadron was first given in Ref. [25] and is given by 

= <^u{x,Pt) + ^L{x,p-r) + ^t{x,Pt) + <^Ll{x,Pt) + ^Lt{x,Pt) + ^Tt{x,Pt), (CI) 

where the contributions ^u{x,Pt), ^Lix,pT) and $7(2;, Pr) are parametrized in the same way as those contributions 
in the correlators that describe an unpolarized, longitudinally or transversely polarized spin 1/2 particle, given in 
Eq. [3] using the notation of Ref. [24] . We update the parametrization for the remaining correlators contributing for 
spin 1 particles, using the same TMD PDFs as in Ref. [2^. This leads for the leading twist TMDs to 



^LLix^pr) 
^lt{x,Pt) 

^tt{x,Pt) 



.flLL{x,pl)S,,,, +ihiLLix,pl) Sr^r.^ > , 



.Pt ■ Sl 



M 2 



JlLT[X,p^) — + giLT{x,PTj<^T '3trM-^75 



M M I 2 

p 

r I 2\l"rati^TT I 2\ uuo PtPt p, 

Jitt{x,Pt) TTn 9iTT[x,PT)et^ 



AP 



TT Vp 



AP 



-75 



'>'1TT\X,Pt)15^i^'^t •^TTpcrJ^+inir^rp(X,PT) — , ^ 



(C2) 



(C3) 



(C4) 
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We note that all polarized quark distributions {g and h functions) in a tensor polarized target, h^j^j^, giLT, h'^LT^ 
^iLT' 5itTj h'f^TT ^^^"^ f^iTT' T-odd. Just as in Ref. [25J, the integrated case is given by 

^LLix) = (C5) 

p 

<^lt{x) = hiLT{xh,e^/S^rpj, (C6) 

$TT(a;) = 0, (C7) 

where 

hiLTix,pl) = h[^j,{x,pl) + hfjy^{x,pl). (C8) 
For identifying the proper rank 1 TMD, it is also useful to define 

hiTTix,pl) = h[j^j.{x,pl) + hf^^^{x,pl). (C9) 



Appendix D: Parametrization of the fragmentation correlator 

The parametrization of the fragmentation correlator for a spin 1 particle, also first used in Ref. [25j, is similar in 
structure to the parametrization of the distribution correlator and is given by 

A(z, fcr) = Au{z, fcr) + Al(z, kj.) + At{z, kr) + All{z, k^) + 1^lt{z, k^) + Arriz, k^), (Dl) 

where the Aij{z,kj,), AL{z,kj,) and AT{z,kj,) are the correlators that describe fragmentation into an unpolarized, 
longitudinally and transversely polarized spin 1/2 particle. For spin only Air{z,kj,) is relevant. The correlators in 
Eq. IDll are with the notation of Ref. [24] at leading twist given by 

Au{z,kr) = S^D,{z,kl)+tHt{z,kl)^'j^, (D2) 
AL(z,k^) = I^G^Uz,kl)SH,75 + HtL(z,kl)SH,^^^, (D3) 



f kr ■ ShT 

Ariz^k^) = GiT{z,kl)^^j^--^j5 + HiTiz,kl)^5$^j 



HMz,k,} M, +^1^^^'^-^ /2' 

ALLiz,kr) = ^D,LLiz,kl)Sh,, + iHtLL{^,kl)Sh,,^^^, (D5) 
ATT{z,k^) = l^D,TT{zX) ^^''^j^f^ +G,TT{zX)^^^^ 



where D^rp^ H^. -fff^^. Gilt, H[]^rp, H^j^j., Gitt, H[j.j. and H^j, are T-odd. The relative sign difference between 
certain corresponding TMD PDF and TMD PFF contributions comes from the definition e"^ — e"'^'"^n+pn_(j, where 
interchanging and n_ gives a relative minus sign [2^. The TMD PFFs Gi(z,A;J), iJi(z,fcJ), Hilt{z, k^) and 
Hitt{z, fcj) are defined in the same way as their TMD PDF counterparts, whereas integrated TMD PFFs are defined 
as 



D 



,{z) ^z'^ j d^k^D„Xz, kl). (D8) 
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